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ABSTRACT 

We present the results of a deep near-infrared imaging survey of the Rosette 
Complex made with the Florida Multi-Object Imaging Spectrograph (FLAMIN- 
GOS) at the 2.1m telescope at Kitt Peak National Observatory. We studied the 
distribution of young embedded sources using a variation of the Nearest Neigh- 
bor Method applied to a carefully selected sample of near-infrared excess (NIRX) 
stars which trace the latest episode of star formation in the complex. Our analy- 
sis confirmed the existence of seven clusters previously detected in the molecular 
cloud, and identified four more clusters across the complex. We determined that 
60% of the young stars in the complex and 86% of the stars within the molecular 
cloud are contained in clusters, implying that the majority of stars in the Rosette 
formed in embedded clusters. Also, half of the young embedded population is 
contained in four clusters that coincide with the central core of the cloud, where 
the main interaction between with the HII region is taking place. We compare the 
sizes, infrared excess fractions and average extinction towards individual clusters 
to investigate their early evolution and expansion. In particular, the size and de- 
gree of central condensation within the clusters appear to be related to the degree 
of infrared excess and mean extinction in a way that suggests that the clusters 
form as compact entities and then quickly expand after formation. We found 
that the average infrared excess fraction of clusters increases as a function of 
distance from NGC 2244, implying a temporal sequence of star formation across 
the complex. This sequence appears to be primordial, possibly resulting from 
the formation and evolution of the molecular cloud and not from the interaction 
with the HII region. Instead, the main influence of the HII region could be to 
enhance or inhibit the underlying pattern of star formation in the cloud. 
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Introduction 



The study of the distribution of embedded stars in Giant Molecular Clouds (GMCs) is 
important for our understanding of the problem of star formation. At the embedded stage, 
these young stars represent the most recent episodes of star forming activity and therefore 
their physical properties should reflect the initial conditions of their parental cores. Near- 
infrared observations of the distribution of embe dded stars in nearby clouds revea l that 
the majority of stars form in embedded clusters (ILada et al.l Il99ll ; ICarpenterl 120001 ) with 
rich clusters (> 100 members) clearly dominating over s mall groups as they contain more 
than 80% of the embedded stellar population in GMCs (IPorras et al.l 120031 ; lLada fc Lada 



20031 . hereafter LL03). However, such surveys have been mostly restricted to nearby clouds 
(d < 500 pc), mainly because the study of more distant complexes is limited by photometric 
depth and resolution of the existing observations. Data from available surveys like 2MASS 
are not sensitive to faint sources (K > 14.3), limiting our ability to sample the young, low 
mass population in distant molecular clouds. This makes it difficult to make meaningful 
comparisons between properties of star forming regions in different parts of the galaxy. 

As part of the NOAO survey program Toward a Complete Near-Infrared Spectroscopic 
and Imaging Survey of Giant Molecular Clouds (P. I.: Elizabeth A. Lada), we chose the 
Rosette Complex in the constellation of Monoceros for a new systematic near-infrared study 
of a molecular cloud outside of the local neighbo rhood. A complete review of p ast and current 
studies on the Rosette Complex can be found in lRoman-Zuniga &: Ladal (120071 ) . The Rosette, 
located in the constellation of Monoceros at a distance of d = 1.6 ± 0.2 kpc (IPark &: Sung 



20021 ; iHensberge et al.ll2000l ). stands out as a very important astrophysical laboratory for the 
study of the early evolution of stars, because of its exquisite layout: at the center of its best 
known feature, the Rosette Nebula, there is a giant OB association, NGC 2244, whose pow- 
erful UV radiation has generated an expanding HII region. It has been suggested that the 
shock front generated by the p hoto-dissociation bub ble could have triggered the formation 
of a family of young clusters (IPhelps fc Ladal |l997J, hereaft er PL97) located in the adja- 



cent, highly structured Rosette Molecular Cloud (RMC; e.g. I Williams et al. 



19951. hereafter 



WBS95), following a process known as Sequential Star Formation (SSF: lElmegreen fc Lada 
19771 ). To test this hypothesis, it is crucial to investigate the influence of the local environ- 
ment (particularly the expanding HII region) on the properties of the young clusters, as it 
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would help us to understand the different initial conditions under which cluster formation 
occurs. 

The Rosette has been previously surveyed in the near infrared. PL97 obtained J,H and 
K photometry with the Simultaneous Quad Infrared Imaging Device (SQIID) at the KPNO- 
2.1m telescope, and were able to identify seven embedded clusters by visual inspection of 
their images. In figure [T] we show the location of th e PL97 clusters^ i n the RMC, in the 
context of the main regions of the cloud identified by lBlitz fc Thadeud (119801 ): The regions 
Al-1 and Al-2, also known as the "Ridge" and the "Central Core" of the cloud, host 4 of 
the PL97 clusters: PL02, PL04, PL05 and PL06. The cluster PL07 is located in region B 
or "Back Core" of the cloud. Finally two more clusters, PL01 and PL03 are located in two 
cores apparently separated from the main cloud, labeled as C and D. All of the clusters 
discovered by PL97 were associated with a luminous IRAS source and a massive molecular 
clump from the study of WBS95, and until now, no other study has been able to determine 
the existence of additional clusters. 



Recently, a spatially uniform catalog from the 2MASS survey was used by iLi fc Smith 



(120051 ) to study the distribution of star forming activity in the Rosette Molecular Cloud. 
They suggested the existence of three general areas of cluster formation, coincident with 
the areas with highes t extinction and temperature of the cloud. The extensions of the star 
formation regions of ILi fc Smith! were determined by eye, roughly tracing the division of 
BT80 and enclosing the main clumps of WBS95 and the locations of the PL97 clusters. 
These general star forming regions could not be compared with the clusters of PL97 because 
the extensions and membership of individual clusters were not determined in the later study. 
Also, due to the low sensi tivity limits of 2MASS and SQIID (K < 14.2 mag), neither the 
study of PL97 nor that of ILi fc Smith! could adequately sample the low mass population at 
the distance of the Rosette Complex and therefore could not determine the fraction of stars 
forming in clusters. Furthermore, the Rosette Complex is located at a low galactic latitude 
(b pd —1.8), and the level of field star contamination in the direction of the cloud is very high, 
making the identification of cloud members even more difficult, specially for low sensitivity 
data. 

In our study we try to overcome these problems by a) using deeper observations of 
the region, capable of detecting sources below the limits previously set by 2MASS and 
b) by using a selected sample of sources with near-infrared exce ss (NIRX) to trace young 
populations with a minimum of non-me mber contamination (e.g . ILi et al.l 119971 ). and c) by 
using the method of Nearest Neighbors (ICasertano fc Hutlll985l . hereafter CH85) to detect 



1 We use from now on the nomenclature "PL01" to "PL07" to refer to these clusters. 



-4- 



surface densities intrinsically larger than the highly crowded field, thus allowing us to detect 
clusters of embedded stars. 

The main goals of our survey are to investigate the distribution of star formation in the 
Rosette, determine the fraction of stars that form in embedded clusters and compare these 
results to those found in other local clouds such as Orion. We also determine the distribution 
of the physical properties (i.p. size and near-infrared excess fraction) of the identified clusters 
in order to study the star forming history and evolution of the clusters. Finally, we aim to 
investigate the distribution of embedded sources across the cloud to investigate the influence 
of the HII region on its star forming properties, and also to test the hypothesis of sequential 
formation. 

The paper is organized as follows: in section [2] we describe our observations and the 
quality of our dataset; in section [3] we explain the methods used to determine the distribution 
of young stars in the Rosette complex, specifically the selection of near-infrared excess sources 
and the use of the nearest neighbor method; in section HJ we present the results from our 
analysis, namely the identification and properties of the embedded clusters; and in section 
M we discuss the implications of our results in the the context of the global properties of the 
Rosette complex along with a summary of our work. 



2. Observations and Data Quality Assessment 

The data were obtained with the Florida Multi Object Imaging Near-Infrared Grism 
Observational Spectrometer FLAMINGOS at the KPNO 2.1m telescope. This project is part 
of the NOAO survey program Toward a Complete Near-Infrared Spectroscopic and Imaging 
Survey of Giant Molecular Clouds. At the 2.1m telescope, FLAMINGOS has a FOV of 
20' x 20', and in the case of the Rosette Complex, a total of 23 FLAMINGOS fields were 
observed during winters of 2001 to 2004 using the J, H and K filters (1.25, 1.65 and 2.2 /im, 
respectively.) The Rosette observations were designed to sur yey the region s of si gnificant 



molecular gas density revealed by the 12 CO and 13 CO maps of iBlitz fc Stark! (119861 . see also 
WBS95) (integrated intensity levels above 20 K-km-s _1 and 0.8 K-km-s -1 , respectively) as 
well as the 25 fim emission map from the IRAS survey (we considered as significant those 
areas with emission above 5.0 MJy-sr -1 ). The result is a very complete coverage of the 
Rosette Nebula and the Rosette Molecular Cloud regions. A map showing the positions 
of the observed fields can be seen in Figure [2j Twenty of the fields are adjacent, while 3 
of them (areas 4, Gl and G2 ) were added to enhance the quality of the observations in 
some particularly interesting regions. In addition, to account for the field contamination, we 
observed 2 near control fields located away from the main molecular cloud emission. These 



- 5 - 



control fields have the same depth as any of our on-source fields, and their locations are also 
indicated in Figure [2j 

All the data of the FLAMINGOS star formation survey were processed with the aid of 



pipel ines built in the standard IRAF Command Language. One pipeline (see iRoman-Zuniga 



20061 ) processes raw data by applying a linearization correction, subtracting dark current, 
dividing by dome flat fields, applying a two pass sky subtraction method, and combining 
reduced frames into image products which are free of geometric distortions, re-sampled to 
twice their size and trimmed into final fr ames of 409 6x4096 pixels. These products were 



then analyzed with a second pipeline (see iLevinei 120071 ). which extracts all possible sources 



from a given field, applies PSF photometry, and calculates an accurate astrometric solution 
and creates final catalogs after matching lists from individual filter observations. Both 
photometry and astrometry are calibrated with respect to 2MASS. 

For the Rosette survey, the completeness limits are K = 17.25, H = 18.00 and J = 18.50 
mag, and are slightly fainter in a few fields observed under optimal weather conditions. The 
completeness limits were determined with artificial star experiments, and apply well within 
75 to 80% of the area covered in each field, but are up to 0.25 mag brighter near the edges of 
the fields, where the photometric sensit ivity of the FLAMING OS detector degraded from the 



center due to a loss in optical quality (lRoman-Zunigall2006l ). This optical degradation also 



affected the shapes of the star images near the edges, and we required to apply a correction 
to the zero point of the photometric calibration. The correction is radial, with center on the 
point of best PSF adjustment (which varied, depending on field and epoch within a radius 
of approximately 100 pixels from position [3000,2400] in the 4096x4096 final image grid), 
and for the Rosette Survey fields it was well defined with a sixth order polynomial form. 
Those stars located at radii where the correction was larger than 0.1 mag were not included 
in the final catalogs. In Figure [3] we show a map of the individual fields and the extension of 
the bad quality areas, which in all cases were smaller than 12%. The zero point correction, 
however, helped to reduce photometric scatter in the accepted areas by an average of 0.02 
mag in each field. Our final mean photometric uncertainties are 0.058 ±0.012, 0.064 ±0.018 
and 0.056 ± 0.014 for J, H and K respectively. 

Figure H] shows the mean photometric uncertainties and their standard deviations per 
magnitude bin as a function of brightness in each band. In this analysis we only used 
catalog stars with photometric uncertainties deviating no more than 3. 0<t from the mean. 
Completeness limits for the analysis were not significantly compromised with these cuts, 
because most of the sources rejected for having large photometric uncertainties were fainter 
than J = 18.50, H = 18.25 and K = 17.50 mag. Additionally, for those sources in which all 
3 bands were required for the analysis, we placed a new limit to the total color uncertainties, 
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(Tj-h = Vcj 2 + &h 2 and oh-k = V^h 2 + <?k 2 , requiring them to be smaller than 0.1 mag. 
This last constraint restricted the photometric errors in each band to be no larger than 0.071 
mag. 

Sources brighter than approximately 11.0 mag suffered of saturation. To account for this, 
we rejected all sources above H = 11.0 from our catalogs (this cut also removed saturated 
stars in bands J and K, which had slightly brighter saturation limits), and replaced them 
with the corresponding 2MASS point sources obtained from the All Sky Release database. 
The total number of sources replaced was 798. After applying the photometric quality cuts 
to the entire database we retained 74,285 sources for analysis. 



3. Analysis 

The Rosette Complex is located at 1.6 kpc from the Sun, more than four times further 
away than regions such as Orion and Perseus, where systematic searches for embedded 
clusters have been previously performed. Furthermore, the Rosette is located at a very 
low galactic latitude (b = —2) and towards the anti-center of the Galaxy (/ = 210), which 
results in a very high density of field sources located in the foreground and background of 
the cloud. These characteristics make it more difficult to detect clusters using traditional 
star count methods. For that reason, our analysis focuses on optimizing the detection of 
young, embedded sources in the cloud. This is done by a) selecting them from a sample of 
infrared excess stars (NIRX) and b) by analyzing their distribution of surface densities using 
a method of nearest neighbors. 



3.1. Infrared Excess Stars 



In the J — H vs. H — K color-color diagram, NIRX stars fall to the right of the band 
defined by the projection of the giant and dwarf sequences along the direction of interstellar 
reddening. This vector of extin ction is defined by a standard reddening law (e.g. ICohen et al. 
198ll ; iRieke fc Lebofskylll985l ). If de-reddened, low m ass NIRX stars would fall back to the 
Classic T Tauri star (CTTS) locus (IMeyer et al.lll997l ). which is defined by large H — K colors 
indicative of dramatic intrinsic extinction due to the presence of circumstellar material. The 
fraction of stars falling i nto this NIRX region has been observed to be fairly large in deeply 



embedded clusters (e.g lLada et al.lll996l ; ICarpenter et al.lll997l ). 



Studies have shown that th e total NIRX fractio n in young clusters decreases with age, as 
a result of inner disk evolution (jHaisch et al.ll2001bl ). For example, the JHK fraction in the 
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Trapezium and NGC 2024 (ages ~1 Myr) is close to 50%, while in slightly older clusters like 
IC 348 (age ~2 Myr) the JHK fraction is around 20%. The deeply embedded population 
of the RMC can be expected to be younger than the exposed OB association NGC 2244, 



which is estimated to ha ve an age of 2 Myr from spectroscopic observations (IPerez et al. 



19891 ; iPark fc Sungj|2002l ). Thus, the JHK excess fraction in the embedded clusters of the 
Rosette can be expected to be between 20 and 50%, significant enough to trace well the most 
recent episode of formation. 

For our study, we define an NIRX star as one with colors that place it 0.1 mag (5 
times the standard deviation of the H — K uncertainty) to the right of the dwarf zero age 
main sequence reddening band (ZAMSRB), and above J — H = 0.47 (H — K) + 0.46, which 
de fines the lowe r limit for the Classic T Tauri Star (CTTS) locus given the listed uncertainties 



of iMeyer et all The first of these two constraints avoids contamination of the sample from 
non-NIRX stars with large H — K uncertainties that could place the stars slightly outside 
the right edge of the ZAMSRB. The second one impedes the inclusion of objects that lie 
in the region located below the CTTS line. This region, just as the one to the left of the 
ZAMSRB, contains detections that are mainly due to high color scatter. Also, i n some cases 



unres olved galaxies (which can pass as stars) with large photometric spreads (ILabbe et al. 



20031 ) or highly inclin ed intrinsic reddening vectors -due to dominant emission from HII 



regions and AO stars (IHeraudeau et al.lll996l )- are located in this region of the diagram 



Combined effects of variable seeing quality over the seasons and variability of the focus 
quality across the wide detector field of FLAMINGOS, resulted in a high dispersion of 
color values at the faint end of our sample, and it cannot be reduced with the zero point 
corrections nor the uncertainty restrictions. The scatter in the H — K colors increases for 
sources fainter than K = 15.75 mag. At this limit, the field object density also increases 
significantly, and the larger the uncertainty in H — K, the more difficult it is to distinguish 
cloud members from background sources. In Figure [5] we show the color-color diagrams for 
stars with K < 15.75 mag and for stars with 15.75 < K < 17.25 mag respectively. In these 
diagrams, the lowest contour level shown represents the mean surface density of objects in 
the color-color space, and each subsequent level represents a step of 1 standard deviation. 
The diagrams show how the scatter for the bright end of the sample is much smaller than 
for the faint end, where we can see a quasi symmetric increase in the dispersion of colors at 
both sides of the ZAMSRB near the most populated areas of the diagram. Our calculations 
indicate that the mean standard deviation of the H — K color dispersion is twice as large 
for the faint end bins (0.11 vs. 0.21 mag, calculated for the entire catalog). The effect is 
also slightly worse for those fields in the survey which were observed under the less favorable 
weather and airmass conditions. 
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Large scatter in photometric colors directly affects the calculation of the number of 
NIRX stars in our survey. We performed a series of Monte Carlo experiments in which 
we simulated the colors of s tars drawn from a model population with an age of 1 Myr 



( ID'Antona fc Mazzitellil 119971 ) located at the distance of the Rosette. We then added color 
errors and extinction to these sources, similar to those observed in the survey fields. We 
found that the number of stars with colors similar to those of NIRX stars due to photometric 
scatter can be up to 5 times larger for stars fainter than K = 15.75. 

Because we are basing our analysis on the detection of infrared excess sources, we limited 
the primary aspect of our analysis, the identification of embedded clusters, to those stars in 
the bright end of the sample (K < 15.75 mag). This sample restriction is rather conservative, 
but assures a more robust detection across the entire survey area. It is important to remark 
that these bright NIRX sources are only helping us to trace the location and rough extension 
of clusters, and although our color uncertainties are high, individual K band magnitudes 
are still good within 0.05 mag (averaged), and thus the sample is still robust enough to 
construct accurate luminosity distributions with a bin resolution of 0.25 mag. Also, the 
photometric depth limit of K — 15.75 implies a sample 1.5 magnitudes deeper than 2MASS 
and is equivalent (for dwarf type stars) to a stellar mass range of 0.09 to 0.18 M & 
for a popula tion of 1 Myr embedded in a cloud with a typical extinction of to 10 visual 
magnitudes (jp'Antona Sz Mazzitelli 1997 ) 



3.2. The Nearest Neighbor Method 



In o rder to optimize cluster identification, we used a variation of the Nearest Neighbors 
Method (ICasertano fc Hutlll985l . hereafter NNM) to determine which objects in our fields 
have surface densities above the crowded background levels. The NNM ha s been shown to 



have success in detecting and outlining the extens ions of embedded clusters f INakajima et al. 



19981 ; iGutermuth et al.ll2005t iFerreira et al.ll2007l hereafter FL07). The generalized form of 
the jth nearest neighbor surface density estimator for a star in a field is simply: 



J 



1 



(1) 



where Dj is the distance from any given star to its jth neighbor. This estimator has 
only one degree of freedom, the fixed number of neighbors, j, used to calculate the local 
density. For a large enough value of j, the fluctuations in the local density estimations due 
to local irregularities will be small, which assures an accurate determination of the extension 
and structure of large systems. However, for the same reason, j also limits the minimum 
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number of grouped stars that can be detected without confusion: if j is much larger than 
the minimum number of stars that define a group or cluster, then small structures could 
be overlooked, as their surface densities will be systematically lower than a defined average 
value. CH85 showed that j = 6 is the minimum number that can assure unbiased detection of 
groups with up to 1000 stars, but is 15% more susceptible to local statistical fluctuations than 
larger values of j. The selection of j for each application will simply depend on the goal: for 
example, in their study of young clusters using 2MASS data, FL07 chose a j = 20 estimator, 



which allowed them to identi 



y clu sters with radii as small as 0.3 pc and a minimum of 20±5 



members. iGutermuth et al.l (120051 ) used j = 5 to determine the membership in cluster fields 



at high resolution in combination with carefully corrected source counts. 

At first, we tried to emulate FL07 by applying a j = 20 estimator to our final RMC 
catalog and our control fields, expecting to be able to detect and delimit the already known 
structures: NGC 2244 and the seven embedded clusters reported by PL97. The result was 
that NGC 2244 could be easily identified as a high density region, along with the clusters 
PL04 and PL05, located in the core of the cloud. The three mid size clusters PL01, PL03, and 
PL07 could be identified too, but their apparent structures were not much more significant 
than some density fluctuations in our fields due to patchy extinction. Finally, stars in the 
smaller clusters PL02 and PL06 had 20th NN densities below the 3a level, and thus these 
clusters could not stand out above the noise levels of the distribution. 

Cluster PL06, which is associated with the B-type proto binary AFGL961, contains large 
quantities of obscuring material near its center, complicating the detection of it s faintest an d 



most embedded members even in carefully constructed near-infrared maps (e.g lAspinlll998l ). 
Cluster PL02 is deeply embedded near the front part of the cloud, where material appears to 
be in contact with the expanding HII region. Likely, for small embedded clusters like PL02 
and PL06, the combination of small number of members, high background contamination 
and non-uniform extinction, hampers the effectiveness of the NNM. Furthermore, in regions 
of high extinction, some embedded clusters could have surface densities intrinsically lower 
than the field so that their members could be undetectable with a large j. Or, random 
patches of low extinction in the cloud could have larger than average field densities and may 
pass as "false" clusters. 

In order to overcome these complications, we applied the NNM only to a JHK color 
selected sample of sources that exhibited a near infrared excess. This allowed us to trace 
preferentially the youngest and most probable members of the cloud. For this reason, the 
number of NIRX sources is expected to greatly diminish outside the molecular cloud, and 
therefore background contamination is reduced, along with false cluster occurrences. Also, 
the number of young sources with NIRX is expected to be larger in deeply embedded popula- 
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tions, and consequently high extinction regions could actually present larger NIRX densities, 
opposite from the "all-star count" case, in which these regions would be overlooked. 

Since it is expected that only a fraction of the members in a given young cluster exhibit 
infrared excess, then our sample will diminish in size, and our j estimator value has to 
be reduced accordingly. We experimented with different values of j, from 6 to 15, and we 
ultimately adopted a value of j = 10 for our analysis. This choice allows us to detect NIRX 
groups of 10 members, which for a typical NIRX fraction of 30 percent (see section 13.11) 
would trace a cluster with ~ 30 members. A preliminary visual inspection of the images of 
the less populated clusters, PL01, PL02 and PL06, indicated that the typical number of stars 
within visible areas of bright nebulosity -usually coincident with embedded cluster cores- 
could be ind eed, around 30. This i s close to the minim um number that defines a cluster 
dynamically (lAdams fc Myersll200ll : lLada fc Ladall2003l ). and below that, any groups could 
at most be considered associations or isolated populations, although they will be harder to 
distinguish from local fluctuations. In principle, a lower value of j would allow us to detect 
those smaller stellar groups, but our tests indicate that analysis with lower values (i.e. j=6 
to 10) do not differ significantly from the case of j = 10, and thus we opted for this, more 
robust choice. 



3.3. Nearest Neighbor Analysis for Infrared Excess Stars 

We detected a total of 1169±34 NIRX sources with K < 15.75. In Figure M we 
show their distributions of 10th Nearest Neighbor distances, D\q and local surface densi- 
ties /iio- The mean value of D w is 3.79 arcmin (1.83 pc) which corresponds to a /iio = 0.2 
stars-arcmin -2 (0.92 stars-pc -2 ) . This limit is also indicated in the plot. 

For the control fields we found 19 sources that had NIRX colors down to a maximum 
brightness of K = 15.75. The mean 10th Nearest Neighbor surface density of the 19 sources 
together is 0.18 stars-arcmin -2 (0.82 stars-pc -2 ), which compares well with the mean value 
0.2 stars-arcmin -2 for the fii distribution in the cloud survey areas. Thus, we considered 
a round value of 0.2 stars-arcmin -2 as a background level for the NIRX population, below 
which we cannot assure that an embedded NIRX source has a density high enough to be 
distinguished from a field source with similar colors. 

The minimum value of the D w distribution in the survey is 0.311 arcmin, which rep- 
resents a density of 29.5 stars-arcmin -2 (136.3 stars-pc -2 ). This value is a good estimation 
of the typical local surface density in the central regions of RMC clusters. The midpoint 
between this minimum distance and the mean is 2.387 arcmin (1.1 pc) which corresponds 
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approximately to yUio=0.6 stars-arcmin 2 (2.77 stars-pc 2 ). It turns out that this value is a 
good estimate of the average embedded cluster size in the RMC. 



4. Results 

4.1. Identification of clusters 

In Figure [7] we show the location of NIRX stars with surface densities higher than 
the background level at 0.2 stars-arcmin -2 . All of the known clusters are traced well by 
these sources, which confirms that they are among the main regions of star formation in 
the Rosette Comple x. In this plot we also show the contours of CO from the study of 



Blitz &: Stark (1986). From now on we will use the nomenclature of the main RMC features 



from lBlitz &: Thadeusl (119801 ). As seen in the figure, most NIRX stars with high NN densities 
are located within the molecular cloud, except for those located in the field of the Rosette 
Nebula. 

In Figure [S] we show a contour level map of the local surface densities calculated with 
the NN method (j = 10). The contours were constructed using a Nyquist sampling box of 
90". Using this map, we define an embedded cluster as a region for which we can observe a 
closed contour at a level of 0.2 stars-arcmin -2 , containing at least 10 NIRX sources. Using 
this definition we found, in addition to NGC 2244 and the seven clusters from PL97, four 
additional areas that arise as significant but have not been previously studied: 

The first new cluster, REFL08 is located in the "Central Core" of the cloud, East of 
cluster PL05 and South of cluster PL04. These two large, neighbor clusters present patchy 
reddening and have a number of members already visible in optical plates (e.g. DSS), thus 
they can be considered to be partially emerged from the cloud. In contrast, REFL08, is 
not visible in any optical plates and has a large number of conspicuously reddened sources. 
Therefore, REFL08 appears to be a third, more deeply embedded substructure in this active 
region of star formation at the center of the cloud, in which clusters PL04 and PL05 are the 
largest clusters. 

The second new cluster, REFL09 is a relatively extended, highly reddened cluster lo- 
cated in the southeastern edge or "Back Core" of the cloud. Along with REFL08, these are 
clear examples of clusters which are located in regions of very high extinction, resulting in 
stellar surface densities comparable or lower than the field if a color selection is not applied. 
However, these clusters contain a large number of heavily reddened sources -easily distin- 
guishable in JHK color composite images- and a significant number of NIRX sources that 
let them stand unequivocally as embedded clusters. 
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The third new cluster is located to the east of NGC 2244, in the region of the cloud 
historically i denti fied as NGC 2237. The existence of this cluster was also suggested by 



Li fc Smith! (120051 ) in their study of 2MASS data. 



The fourth cluster, REFL10, is much less prominent. It is located North of NGC 2244, 
and despite of its low NIRX surface density compared other star forming regions in the 
complex, its 13 NIRX sources make it qualify as a cluster. REFL10 is an example of cluster 
that would not be significant if the NN distribution of NIRX sources was analyzed with a 
large value of j, or if a color selection was not applied. 



4.2. Properties of Clusters 

Once the clusters were identified, we proceeded to investigate their basic properties 
and how are they related in the context of the evolution of the star forming complex. For 
each cluster we determined sizes, centers, number of members and infrared excess fractions. 
These cluster properties were determined individually, isolating regions centered at each 
cluster density peak. The area of these individual cluster analysis regions varied from 25 to 
120 arcmin 2 , depending on the apparent extension of the cluster. This way we were able to 
observe in detail structures occupying areas between 6 and 60 arcmin 2 . 



4-2.1. Determination of cluster properties 

Cluster boundaries are defined by the 0.2 stars-arcmin -2 contour level in our NIRX NN 
density plots. All sources with K < 17.25 mag within this contour level are considered 
potential cluster members. Cluster sizes are determined by measuring the area, Ap, within 
the cluster boundary and calculating the equivalent radius, R eq = a/ Ap/ir. 

The density peak, X d , and the core radius, R cor e, of each cluster were calculated with 
the formulation of CH85. The density peak defines the natural center of a cluster, and is 
simply the density weighted average of the star positions in a given region: 



EiW) 

Similarly, the core radius, R cor e is defined as the density weighted average of the distance 
of each star to X d j\ 
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In Table [U we present the center coordinates (density peaks), core radii, equivalent 
radii, and average extinction values for each cluster. The latter were calculated from the 
H — K colors. We also list the number of NIRX sources to K < 15.75 mag and the 
corresponding excess fraction. The excess fraction was calculated by comparing the number 
of IR excess sources to the number of cluster members. The number of cluster members 
was determined by subtracting the expected number of background field stars from the 
number of stars having K < 15.75 within the cluster boundary. The number of field stars 
was determined for each cluster by scaling our off fields to the cluster areas and applying 
appropriate extinction corrections. 

In Figures M to [2D] we show images of each of these cluster regions, as well as the color- 
magnitude diagrams, color-color diagrams and stellar radial profiles for all stars within the 
cluster boundaries. In the images, we marked NIRX sources brighter or equal than K = 15.75 
mag with 'X' symbols. In the K vs. H — K color-magnitude diagrams we plot, for each 
cluster, all of the stars inside the corresponding 0.2 stars- arcmin -2 contour, and we marked 
with 'X' symbols those objects with infrared excess. We also plot the ZAMS locus and a 
PMS evolution isochrone of 1 Myr, both shifted to the distance of the Rosette. The indicated 
extinction vectors corresponding to 3 times the mean value (Ay) in the cluster analysis box. 
Stars falling to the right of the 1 Myr isochrone are clearly affected by extinction towards 
the line of sight of the cluster, revealing their highly embedded nature. In the J — H vs. 
H — K color-color diagrams, the same stars are located above the dwarf and giant star 
sequences along the reddening bands, with the NIRX sources located to the right of the 
ZAMS reddening strip. Those objects located at or near the zero age sequences are most 
probably foreground stars or evolved cloud members. 

The radial density profiles in the fourth panels of each cluster fi gure were constructe d 



using annuli of decreasing width, each enclosing the same area (see e.g. iMuench et al.l (120031 )) 
and centered on the weighted density peaks defined above. For these radial profiles we took 
into account all stars down to K — 17.25 inside the analysis boxes. In the plots we indicate 
the equivalent and core radii defined above. Clusters PL02, PL06 and REFL10, have the 
lowest surface densities and their profiles are poorly defined, but the rest present well defined 
radial profiles, which show well extended tails. In four cases — clusters PL01, PL04, PL07, 
REFL09 — the profiles show secondary peaks suggestive of structure. In the case of the 
Nebula clusters NGC 2244 and NGC 2237, their radial distribution profiles decline slowly, 
which suggest a diluted core peak and a rather extended structure. 
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4-2.2. Comparison of cluster properties 

Under the assumption that the NIRX sources trace well the extent of the clusters, the 
core and equivalent radii are good estimations of the total sizes of the clusters. The core 
radii, R cor e of the Rosette clusters have a range of 0.69 to 4.17 arcmin (0.32 to 19.39 pc), 
with an average of 2.0±1.0 arcmin (0.93±0.47 pc). Their equivalent radii, R eq , vary from 
1.6 to 4.9 arcmin (0.74 to 2.28 pc), with an average of 3.1±1.0 arcmin (1.44±0.47 pc). The 
distributions of these size estimates are shown in the top panel of Figure [2TJ 

In the bottom panel of Figure |2~T1 we show the distribution of the ratio r = R core /R eq , 
which has a nominal peak at r = 0.7. In the study of FL07, where NN analysis was applied 
to a large sample of known cluster fields, they found that embedded clusters separated in 
two groups by their value of r: those with ratios below 0.5 are denominated as centrally 
condensed (C-type clusters), while those with values closer to 1.0 are denominated as having 
a flat profile (F-type clusters). We found that only two of the RMC clusters, PL01 and PL03, 
have C-type ratios, and for three of the clusters, PL06, REFL10 and NGC 2237, cores could 
not be easily determined and, for purposes of comparison, were given a ratio of r = 1.0. 

The distribution of r ratios is rather different from the one observed by FL07 in their 
larger and more heterogeneous sample of clusters: they found a larger fraction of compact, 
embedded clusters (44 percent of the sample) with well defined cores, while in our sample, 
centrally condensed clusters appear to be scarce. However, we note that the NN analysis of 
FL07 was not applied on color selected samples, but on all sources detected in a given field. 
From the observed distribution of cluster sizes and core-to-total size ratios it appears that 
the Rosette clusters have preferentially extended profiles. 

As listed in Table [IJ the observed NIRX fractions of the nebula clusters, NGC 2244 and 
NGC 2237 are 12 ± 3 and 15 ± 2%, respectively, and therefore they are significantly smaller 
than those of the embedded clusters, which vary from 18 ± 3% in cluster PL05 to 76 ± 9% in 
cluster REFL09. Also, the average NIRX fraction in the molecular cloud areas is 41 ± 6%, 
i.e almost three times as large than in NGC 2244 and NGC 2237. Curiously, the fraction 
observed for REFL10, which is also located at the nebula is much closer to the average of 
the embedded clusters and actually slightly larger than those of clusters PL01, PL04 and 
PL05. 

In Figure [22] we show the relation between the NIRX fraction, the median extinction 
value and the equivalent radii of clusters. Extinction and cluster size are negatively correlated 
with a Pearson coefficient r\ = —0.55, while extinction and NIRX fraction are positively 
correlated with a Pearson coefficient r] = 0.57. Consequently, as shown in Figure [231 cluster 
equivalent sizes and core sizes are both negatively correlated with NIRXF, with Pearson 
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coefficients r\ = —0.61 and r\ = —0.54 respectively. 

We also found that the ratio r was negatively correlated with NIRXF for clusters with 
a distinguishable core, this correlation is weaker (rj = —0.48) but is similar to the relation 
found by FL07 for a sample of nearby clusters. The significance of these correlations in the 
context of the structure of evolution of the RMC clusters is discussed in section |5j 



4.3. The Fraction of Stars in Clusters 

Since NIRX sources trace the youngest population of stars within a star forming region, 
we used the distribution of NIRX sources to investigate the fraction of stars located and 
forming in clusters. The fraction of NIRX sources in clusters was calculated by simply 
counting the number of NIRX sources within the cluster boundaries and comparing this 
value with the number of NIRX sources outside the cluster boundaries, correcting both for 
background contamination. We first considered the entire area of our survey, which provided 
us with an estimate of the number of stars that are currently located within clusters across 
the entire complex. Since the most recent episodes of star formation should be associated 
with the molecular gas, we next considered only the area within the molecular cloud, which 
provided us with an estimate of the fraction of stars currently forming in clusters. 

For the first determination, we calculated that the total area covered by the FLAMIN- 
GOS survey is 7308 sq. arcmin (1579 sq. pc) while the 11 clusters identified with our analysis 
occupy a total surface of 390 sq. arcmin (84 sq. pc), i.e. 5.3 percent of the total survey area. 
The total number of NIRX sources detected is 1169 ± 34, out of which 549 ± 23 are located 
within the boundaries of the clusters, and 620 ± 25 are outside the boundaries. In order 
to correct for the background contamination, we scaled the the number of NIRX sources 
observed in the control fields by a factor equal to the area of the molecular cloud minus the 
areas of the clusters and then divided by the area of the control fields. This expected number 
of NIRX sources due to background contamination was estimated to be 261 ± 17. The back- 
ground correction for individual clusters is always smaller than the Poisson uncertainties but 
it was also applied. Using these corrections, the number of NIRX sources outside of clusters 
was 359 ± 19, and inside clusters 531 ± 23. Therefore, the fraction of NIRX sources located 
in young clusters across the Rosette Complex is 60 ± 5%. 

The area of the molec ular cloud was defi ned by the 0.8 K-km-s -1 integrated intensity 
contours of 13 CO emission Jrlever et alil2006h. corres ponding to Ay = 3.5 ± 0.6 mag using a 



recently derived conversion ratio (IPineda et al.l 120071 ). This area equals 2684 sq arcmin (580 



sq. pc.) The 9 clusters located in the cloud occupy a total of 242 sq. arcmin (52.3 sq. pc) or 



-16- 



about 9% of the RMC The total number of NIRX sources detected in the molecular cloud 
is 589 ± 24, out of which 436 ± 21 are located within the cluster boundaries and 153 ± 12 
are outside the boundaries. After correcting for field contamination, the number of NIRX 
sources located in clusters is 429 ±21 and outside clusters is 73 ± 8. Therefore, the fraction 
of NIRX sources currently forming in embedded clusters in the Rosette Molecular Cloud is 
86±4%. 

Interestingly, a total of 208±15 NIRX sources are contained in clusters PL04, PL05, 
PL06 and REFL08, all located at the "Central Core" of the cloud. This corresponds to 
48±3% of the total number of embedded sources. Therefore, approximately half of the 
recent births in the RMC have occurred at the most dense region of the cloud. 

4.4. Distribution of Sources with respect to the Rosette Nebula 

We analyzed the spatial distribution of those NIRX sources with 10th NN densities 
higher than the mean as a function of the distance to the center of NGC 2244. To do this, 
we determined the surface density of sources within concentric rings centered on NGC2244 
having a width of 1.0 pc. For the first 11 pc we were able to use complete rings, but further 
away, we had to limit the angular extent of the rings to conform with the irregular shape of 
the survey boundaries. The surface density of sources in each ring or segment was calculated 
by dividing the number of sources in each segment by its area and then diving by the surface 
density inside the first parsec circle. 

This normalized radial NIRX source density distribution is shown in Figure [24] (top 
panel). In the plot we indicated the approximate locations of clusters and main cores of 
the molecular cloud, measured in projected distance from the center of the Nebula. We also 
indicate the location of the main regions of the RMC, the "Ridge", the "Central Core", and 
the "Back Core" with respect to the Rosette Nebula. 

In the bottom panel of the figure, we indicate the average fraction of NIRX sources 
in clusters as a function of distance, averaged in five groups of clusters which are roughly 
traced by the distribution shown in the top panel. The cluster groups used are as follows: 
first group includes NGC 2244, NGC 2237 and REFL10; second group includes PL01 and 
PL02; third group includes PL04, PL05 and REFL08; fourth group includes PL03 and PL06; 
finally, group 5 includes clusters PL07 and REFL09. 
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4.5. Low density population 

In order to investigate the low density star forming component in the cloud, we closely 
examined the region between the central core and the back core. This region corresponds 
to Field 09 of our survey. This field was chosen because the seeing and observing conditions 
for it were particularly good. Specifically, the average scatter of colors down to K — 17.25 
remains below 0.109 mag (similar to the average for bright end bins) across the whole field 
(see Figure 125]) . This superior quality was achieved because the southeastern quadrant of 
the field, which in other fields presents high stellar profile distortions, overlaps with the 
good quality northwestern quadrant of field Gl. Selective averaging of photometry in these 
overlapping quadrants, permitted the reduction of the total color scatter. 

No clusters were found in this field, but a small group of NIRX sources coincide well 
with a roughly filamentary core in the cloud located southwest of the PL06 cluster region (see 
Figure [7]). Our NNM analysis does not produce a new cluster identification for this core, but 
we counted the number of NIRX sources in the whole field down to brightness limits of 15.75, 
16.25, 16.75 and 17.25 mag. We counted the number of NIRX in the control fields at these 
same magnitude limits. In order to have make a conservative estimate, we limited the counts 
to a circular area with a radius of 1500 pixels centered on pixel position [3000,2400]. This 
delimits the area of best photometric quality after the polynomial correction, as described 
in section [21 

We compared these NIRX counts in field 09 with the expected number of field NIRX 
sources at these same limits in an equivalent area of high photometric confidence in the 
control fields, adding again, a mean extinction level of 5.0 mag to account for stars that 
would have excess colors after reddening. A second comparison was done by averaging the 
number of observed NIRX sources in equivalent photometric confidence circular areas of 
survey fields 03, 13, 14 and 15. These fields are located in regions of the complex where 
the CO emission from the molecular cloud is lower than average, and thus they can help to 
determine if the NIRX sources of field 09 could indeed be a local enhancement near the core 
of the RMC. 

In figure [261 we compare the counts described above in the form of cumulative his- 
tograms: apparently, field 09 surpasses significantly the expected number of NIRX from the 
control field and the off-cloud areas, with counts 4.8, 5.6, 3.5 and 2.2 times larger than those 
in the control field field at K < 15.75, 16.25, 16.75 and 17.25 respectively. The counts are 
also higher than the averaged counts of survey fields 3, 13, 14 and 15 by factors of 1.3, 1.5, 
1.6 and 1.8 respectively. This suggests that the region of the cloud observed in field 09 could 
have a significant number of young source candidates not associated with clusters. Further- 
more, JHK color images of selected regions of the field show a number of highly reddened 
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sources. Some of these sources have thin nebulosities, and coincide or are located close to 
stars with near infrared excess emission. These sources are not included in our infrared 
excess catalog because they were not detected in our J band observations and therefore we 
do not have three colors to determine if they exhibit excess emission. 

In addition, we counted the number of sources with no J band photometry and H — K > 
1.5 mag at the same K brightness limits mentioned above. It is not possible to assure that 
these sources belong to the cloud, especially for K > 15.75, where the scatter increases and 
the contamination by background galaxies worsens. However, we noticed that the number of 
these sources which locate outside the main CO emission regions of the molecular cloud was 
very small; the majority of them were in fact located within the cluster areas (see Figure 
1271) . Therefore, these red sources could be tracing the location of young sources too. Also, 
no stars with H — K colors this large were found in the control field. In Figure [23 we show 
the distribution of the sources with H — K > 1.5 in the survey areas down to K < 17.25 
mag respectively. In table [2] we show the number counts of H — K > 1.5 sources, separating 
those located outside the 20.0 K-km-s _1 CO contours, (Ay ~ 3.5 mag) and then dividing 
the ones located inside the contours into those within and outside the cluster areas. After 
correcting by background (in this case by subtracting the scaled number of in-survey /off- 
cloud sources) we estimate that the average fraction of objects with H — K > 1.5 located 
inside the molecular cloud projected area but outside the clusters is approximately 21±4%. 
The fraction is slightly larger than the one obtained from counts of NIRX sources (14±4%), 
possibly due to the fact that our red source counts are less restricted, as we cannot limit 
their range in J — H. Counts of sources with large H — K values have been us e d before to 



ident ify highly embedded clusters with good results (e.g. iHomeier fc Alvesll2005l ; lOjha et al. 



20041 ); in our case, our maps show that large H-K color sources trace well the regions we 
identified as clusters, but also show a population at the 10 to 20 percent level that might be 
distributed along the rest of the cloud. 



5. Discussion 

5.1. Clusters as a Dominant Mode of Formation 

Observations of nearby molecular clouds such as Orion, Perseus and Monoceros suggest 
that most stars in these clouds form in embedded clusters (LA91, CA00). In the RMC, 
approximately 86 percent of the present day star formation occurs in embedded clusters, 
in excellent agreement with the previous results for nearby clouds. This indicates that 
cluster formation is also the dominant mode of star formation in clouds beyond 1 kpc. 
Cluster formation appears to be well distributed across the Rosette Complex but the clusters 
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themselves are confined to a rather small fraction (5 percent) of the total area of the cloud. 
Comparison of the distribution of embedded clusters with the distribution of molecular gas 
(WBS95), reveals that the embedded clusters are associated with the most massive molecular 
clumps in the cloud, which are also the mos t dynamically evolved. This is similar to what 
was found in the Orion B cloud (jLadal 119921 ). 

While the most massive clumps in the Rosette are responsible for the production of the 
embedded clusters, then some of less massive clumps could be associated with the distributed 
population. The number of sources exhibiting near-IR excess in field 09 and the number of 
sources with H — K > 1.5 across the cloud could indicate that a distributed population 
of young stars formed in addition to the cluster population. However, the formation of 
this distributed population could account for no more than 20 percent of the total star 
formation within the RMC Characterizing the distributed population is difficult due to 
several observational constraints. For example, the large photometric scatter of faint sources 
along with the field star contamination restrict our study to the brighter (likely higher mass) 
sources, which are expected to be less numerous than the fainter (likely lower mass) ones. 
Spectroscopic observations and mid-infrared imaging would help to establish the number 
of highly reddened sources outside the cluster regions. Also, the nature the distributed 
population in molecular clouds is somewhat unclear as discussed by CA00, precisely because 
these sources permeate through large areas of the molecular cloud in and between regions of 
cluster formation and thus do not have a clear origin. Such a population could be the residual 
of a slightly older generation of star formation, a sub-product of cluster formation resulting 
from mass segregation or ejection or belong to a population located in the foreground or 
background of the cloud. 



5.2. Rapidly Evolving Clusters 

The mortality rate of embedded clusters in the Milky Way is quite high. Less than 
10 percent of clusters survive their emergence from molecular clouds and live longer than 
lOMyr and less than 4-7 percent survive to become bound clusters the age of the Pleiades 
(LL03). This early disruption is most likely linked to low star formation efficiency and the 
rapid destructio n of the molecular clouds. Similar results are found in other galaxies (e.g. 



Fall et al.ll2005h . 



The observed distribution and properties of the embedded clusters in the Rosette Com- 
plex may be revealing some of the first signs of such early cluster evolution. The clusters 
in the Rosette appear to be, on average, slightly larger than those located in other clouds: 
the cluster equivalent radii range from ~ 0.75 to 2.3 pc, with an mean of 1.44 and a median 
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of 1.46 pc. The largest cluster in the complex is NGC 2244, with an equivalent radius of 
2.3 pc. For comparison, the sizes of clusters listed in the catalog of LL03 range from ~ 0.3 
to 3.8 pc, but have a mean of 0.8 pc and a median of 0.62 pc. Only a few clusters in this 
catalog, such as MonR2, Gem4, NGC 2282 and Trapezium/ONC, have radii larger than 1.5 
pc. 

It is unlikely that the Rosette is systematically forming larger than average clusters. 
Instead, the large cluster sizes may be due to dynamical evolution of the cloud leading to the 
beginning of cluster expansion and eventual dispersal. Indeed our data suggest that cluster 
sizes in the Rosette may be related to cluster evolution in several ways. First, we found 
a negative correlation between extinction and cluster sizes (Fig. [22]) . indicating that the 
most extended clusters are the least embedded and are the ones that are emerging from the 
parental molecular material. Cluster sizes are also inversely proportional to their infrared 
excess fraction (Fig. [231 ) It has been shown that the infrared excess fraction in young 
clusters is initially high and decreases rapidly as a function of c luster age such that clusters 



with ages > 5 Myr have no or extremely small excess fractions (lHaisch et al.ll2001bl ). In the 
Rosette, excess fraction decreases with increasing cluster size, implying that clusters expand 
relatively quickly (in a time scale comparable to the T Tauri phase). Furthermore, analysis 
of the density structure of the embedded clusters may also indicate the beginnings of cluster 
evolution. For example, even though some of the Rosette clusters exhibit a well defined core 
from our NN analysis, others appear to be extended well beyond this core. Indeed, the ratio 
of the core size to the total size, r, is 0.5 or above for 6 out of 8 young clusters with defined 
cores (Fig. |2"T1) indicating that most of the Rosette clusters do not have prominent central 
condensations. In addition, r decreases with increasing infrared excess fraction (see Fig. 123]) . 
A similar correlation was found by FL07 in their analysis of local clusters. In the Rosette, 
clusters with r > 0.5 have an average infrared excess fraction of 27 percent whereas clusters 
with r < 0.5 have an average IRX fraction of 52%. These results suggest that clusters form 
as compact units and later become less centrally condensed. This initial expansion appears 
to start quickly after formation and likely develops within a few Myr, favoring a scenario of 
rapid disruption and cluster dispersal. 



5.3. Influence of the HII Region and Sequential Star Formation 



The layout of the Rosette Complex led WBS95 and PL97 to suggest that the formation 
of the clusters was triggered by the expansion of the nebula into the cloud , resembling the 
picture of sequential star formation proposed by lElmegreen fc Ladal (Il977l ). In this model, 
the shock front from an expanding HII region would trigger the formation of a new OB 
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association in an adjacent molecular cloud which would later trigger a third group and so 
on. We can use the distribution of young stars and clusters in our survey to investigate what 
effect, the HII region has had on the star forming properties of the cloud. 

Star formation in the Rosette appears to be concentrated into 4 main areas of the 
complex as illustrated in Figure 24. The most prominent is the Rosette Nebula itself, which 
contains the clusters NGC 2244, NGC 2237 and REFL10. The next concentration of star 
formation is located between 10 to 20 pc from the center of NGC 2244, corresponding to the 
RMC ridge. This region is located the closest to and overlaps with the ionization front of 
the HII region and contains clusters PL01 and PL02. It has the lowest levels of star forming 
activity. Star formation may be inhibited in this area due to the effects of the ionizing 
radiation. Hot ionized gas might accelerate the evaporation of gaseous material surrounding 
a forming cluster, resulting in the end of the star forming process. 

The third concentration of young stars is located 20 to 30 pc from NGC 2244 and 
contains clusters PL04, PL05, REFL08 and PL06 at the RMC "central core" and PL03 
located in a gas clump (core D) outside the main area of the cloud. The four clusters in 
the central core account for 48% of the present day star formation in the molecular cloud 
while only covering roughly 4% of the cloud area. These cores are associated with four 
of the most massive cores in the cloud (WBS95). Similarly high concentrations of star 
formation have been observed in other molecular clouds such as Orion, W3/W4/W5 and 
Perseus, where approxima t ely half of the young stars are found in only 1 to a few embedded 
clusters (ILada et al.lll99ll ; ICarpenter et al.ll2000l ; IjQrgensen et al.l 120061 ). The RMC central 
core appears to coincide with the edge of the Rose tte Nebula as shown by i mages of radio 
cont inuum and 2 5 um IR AS emission. The studies of lCelnikl (119851 ) ; ICox et al.l (11 9901 ). WBS95 
and iHeyer et al.l (120061 ) agree that the nebula and molecular cloud have their most intense 
interaction at this region and this may indicate the recent passage of the shock front from 
the nebula. The interaction with the shock front may have triggered or stimulated the 
production of stars here resulting in the second most active region of star formation in the 
complex. 

The last concentration of star forming activity is at the "back core" of the cloud. This 
region contains clusters PL07 and REFL09. These two clusters appear t o have formed in a 
region of the cloud located beyond the main interaction with the nebula (IHeyer et al.ll2006l ). 
Given their distance from the shock front, it is unlikely that the formation of these clusters 
was triggered. This reinforces the suggestion by PL97 that star formation in this part of the 
cloud occurred spontaneously. 



Surprisingly, we do find the first evidence for a possible temporal sequence in star forming 
events in the Rosette. The average cluster infrared excess fraction increases as a function of 
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distance from NGC 2244, suggesting that the embedded clusters are progressively younger 
the further they are from the Nebula. The existence of this apparent age sequence indicates 
that star formati on in the Rosette did indeed take place sequentially in time, but not in the 
way proposed by lElmegreen &: Ladal . Instead, it appears that the overall age sequence of 
cluster formation is independent of any interaction with the expanding HII region but rather 
may be primordial, possibly resulting from the formation and evolution of the molecular 
cloud itself. While the existence of the HII region cannot be responsible for the sequence 
of cluster ages, it does appear to have a significant impact on the underlying sequence by 
either enhancing or inhibiting the star forming process. 



5.4. Summary 

1. In this paper we present results of a deep near-IR survey of the Rosette Complex. The 
survey was made with the wide field imager FLAMINGOS at the Kitt Peak 2.1m telescope 
and covers all the main areas of the Rosette Nebula and the Rosette Molecular Cloud 

2. We analyze the distribution of young stellar sources in the complex by estimating 
the surface densities of objects with infrared excess, using a variation of the nearest neighbor 
method, which allows us to find clusters with minimum populations of ~ 30 members. 

3. We confirmed the existence of the seven embedded clusters found by visual inspection 
by PL97 and found four more clusters in the complex. Two of these clusters are deeply 
embedded in the molecular cloud and the other two are located west of NGC 2244 in the 
Rosette Nebula. 

4. The young cluster population accounts for 60 percent of the stars in the complex, 
and approximately 86 percent of the youngest stars, located in the molecular cloud. This 
implies that the majority of stars in the Rosette form in embedded clusters, similar to nearby 
clouds. . 

5. The sizes of clusters in the Rosette Complex appear to be anti-correlated with their 
mean extinctions and infrared excess fractions, which suggests that clusters form as compact 
units then expand shortly after formation. The timescale for this process is similar or even 
shorter than the T Tauri phase, as evidenced by the significant NIRX fractions and deeply 
embedded status. 

6. The distribution of young clusters suggests a division of star formation into four main 
regions coincident with the largest features of the cloud: the first is the Rosette Nebula, which 
contains the oldest clusters, the second is the ridge, located near to the ionization front of the 
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HII region, the third is the central core of the molecular cloud, where the main interaction 
between the nebula and the cloud is located, and approximately 50 percent of the young 
stellar population was produced. The fourth region of formation is at the back core of the 
cloud, located furthest from the Nebula. 

7. The averaged infrared excess fraction of these regions appears to increase as a function 
of distance from the Rosette Nebula, which is suggestive of a sequence of cluster ages across 
the cloud. Rather than being triggered by the HII region, this sequence appears to be 
primordial, possibly resulting from the formation and evolution of the molecular cloud. The 
HII region appears to enhance or inhibit the underlying pattern of star formation in the 
cloud. 
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Fig. 1. — Integrated intensity 13 CO map of the Rosette Molecular Cloud from the survey of 
Heyer et al. (2006). The locations of embedded clusters identified by Phelps & Lada (1997) 
are indicated by solid dots. The principal regions of the Rosette Molecular Cloud, identified 
by Blitz & Thaddeus (1980), are labeled. 
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ROSETTE COMPLEX 




HJ010 WviJW *7^MU JtWJ «vO»i 

FLAMINGOS Gkinr Molecular Cloud Survey 

Fig. 2.— A map of the FLAMINGOS-NOAO-University of Florida Rosette Molecular Cloud 
Survey. The boxes delimit individual FLAMINGOS fields (20 x 20' after trimming) over a 
halftone image (in logaritmic scale) of the IRAS 25/xm emission in the region. The labels at 
the left side of each box refer hereafter to the fields detailed in table 1 and the text. Light 
solid contours represent the extension of the Rosette Molecular Cloud in CO integrated 
emission from the survey of Blitz & Stark (1986) in steps of 20.0 K-km-s -1 . Crosses mark 
the centers of known embedded clusters from the previous study of Phelps & Lada (1997). 
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Fig. 3. — The extension of the areas of acceptable optical distortion are marked for each 
field as circles with radii equal to the center of the maximum bin at which the polynomial 
correction to the zero points (see text and figure 5) can be applied within the detector. This 
effect varies by field (size of the acceptable area) and filter: the solid, dotted and dashed 
linestyle circles represent the tolerance radii for J, H and K respectively. 
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Fig. 4. — Intrinsic quality of the Rosette Molecular Cloud survey. The points mark the 
mean value of the photometric error in a given magnitude bin, with the size of one standard 
deviation indicated by the error bars. The horizontal dashed lines marks the zero level and 
the sensitivity level, estimated at 0.1 mag. Vertical dotted lines indicate the sensitivity 
limits, marked as the bin at which the fiducial curve generated by the mean values crosses 
the 0.1 line. The dash-dotted line represents the 3-a level of the fiducial errors; any star in 
our catalog with errors higher than this levels were rejected from the analysis. 



Fig. 5. — Contour level color-color diagrams for stars in the FLAMINGOS RMC survey 
divided in two ample groups of brightness. Both diagrams were constructed with a Nyquist 
box size of 0.1 mag. Left: distribution of colors for 'bright' stars within 5.0 < K < 15.75 mag, 
Right: distribution of colors for 'faint' stars within 15.75 < K < 17.25 mag. The solid line 
contours start at the mean level (360 and 466 dex~ 2 for the left and right panel, respec- 
tively) with subsequent steps of 1 sigma (1770 and 1910 dex~ 2 for the left and right panel, 
respectively). 
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Fig. 6. — Nearest Neighbor distributions for bright NIRX stars. The top panel shows the 
distribution of 10th neighbor distances. The bottom panel is the distribution of 10th neighbor 
densities. In the top panel line A indicate the limit of distances shorter than 1.0 pc, while the 
dashed line B indicates the midpoint value at 1.83 pc. At the bottom panel the equivalent 
limits in density space are also indicated. 
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Fig. 7. — The location of NIRX stars in the Rosette survey with brightness K < 
15.75 mag. The "plus" symbols are NIRX stars with 10th neighbor densities higher than 
0.2 stars-arcmin -2 , while black dots are stars with densities below 0.2 stars-arcmin -2 . We 
also indicate the expected centers for the known clusters NGC 2244 and PL01-PL07. The 
contours indicate the baseline level of 13 CO emission of 0.8 K-km-s _1 that we used to define 
the extension of the main molecular cloud regions. The solid thin line indicates the limits of 
the survey coverage. 
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Fig. 8. — Final identification of clusters in the Rosette Complex. The colored scale contours 
indicate levels of surface density from the 10th NN analysis. The contours were constructed 
using a Nyquist box sampling of 1.5 arcmin. The dot line contours indicate the baseline level 
of 13 CO emission of 0.8 K-km-s _1 that we used to define the extension of the main molecular 
cloud regions. The solid thin line indicates the limits of the survey coverage. 




Fig. 9. — Top Left: K band image (white symbols indicate locations of NIRX sources), Top 
Right: control magnitude diagram, Bottom Left: color color diagram and Bottom Right: 
Radial Density Profile for the area corresponding to cluster PL01. See text for explanation. 
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Fig. 10- 



Same as Figure [9j for cluster PL02. 
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Fig. 11- 



Same as Figure [9j for cluster PL03. 
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Fig. 12- 



Same as Figure [9j for cluster PL04. 
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Fig. 13.- 



Same as Figure [9j for cluster PL05. 
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Fig. 14. — Same as Figure [9j for cluster PL06. 
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Fig. 15.- 



Same as Figure [9j for cluster PL07. 
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Fig. 16.- 



Same as Figure [9j for cluster REFL08. 
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Fig. 17.- 



Same as figO for clusters REFL09. 
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Fig. 18. — Same as Figure [9j for cluster REFL10. 




Fig. 19.— Same as Figure El for cluster NGC 2237. 
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Fig. 20- 



Same as Figure El for cluster NGC 2244. 
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Fig. 21. — From top to bottom: distribution of core radii, equivalent radii and core to 
equivalent radii ratios for the Rosette clusters 
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IRXF 



Fig. 22. — Top: cluster equivalent radii vs. median extinction. The solid line and dotted 
lines represent a least squares linear fitting ±1 rms deviation. Bottom: NIRX fraction 
(K < 15.75 mag) vs. median extinction 
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Fig. 23. — Clusters equivalent sizes R eq , cluster core sizes R CO re, and core to total size ratios, 
r as a function of NIRX fraction (top, center and bottom panels respectively). The solid 
line and dotted lines represent a least squares linear fitting ±1 rms deviation. 
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Fig. 24. — Top panel: Distribution of NIRX sources density a function of distance from 
the center of the Rosette Nebula (NGC 2244). The counts are made in sectors of 1.0 pc in 
length and counts in each sector have been scaled and normalized to the area and counts in 
the central 1.0 pc circle in NGC 2244. Labels indicate the approximate locations of clusters 
described in this paper, as well as the main 'regions' of the complex, defined by Blitz & 
Thadeus (1986). Bottom panel: averaged NIRX fractions in each of the cluster groups 
defined from the top plot appear to increase with distance from the Rosette Nebula. 
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Fig. 25. — Distribution of color errors in Field 09 of the survey. The color scatter is low 
in this area due to the combination of bad quality areas with the good quality areas of 
overlapping observations. The dashed horizontal lines mark the 0.0 and 0.109 mag limits. 
The dashed vertical line marks the completeness limit of the observations at K — 17.25 mag. 
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Fig. 26. — Cumulative Counts of NIRX sources in field 09 of the survey. The solid white 
bar histogram indicates the number of NIRX sources in the field down to brightness limits 
of K=15.75,16.25,16.75 and 17.25 mag inside the area of best PSF photometry quality. The 
dotted line bar histogram indicates the average of NIRX counts inside the same circular area 
for fields 03, 13, 14 and 15, all located in areas mostly devoid of strong molecular hydrogen 
emission. The shaded histogram indicates the scaled number of "expected" background NIRX 
sources, averaged from counts in the control fields (also in the best PSF circular areas), and 
adding an uniform extinction value of 5.0 mag. 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



IRX counts F9 




Rosette Complex, J = null, H-K>1 .5, K<17,25 
5.60000 I 



'n 




33.2000 38.BO00 38.4000 38.0000 37.6000 37.2000 

Right Acens'ofT ( J 2000) 

FLAMINGOS GMC Survey. HOAO/UF 



Fig. 27. — Distribution of sources with no detection in J and H — K > 1.5 (circled dot 
symbols), limited to brightness of K = 17.25. The thin solid line contours indicate 13 CO 
emission in steps of 1.0 K-km-s -1 . 
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Table 1. Young Clusters in the Rosette Complex 



Cluster 
ID 


RA 


DEC 


Rcore 

[P' 


Requiv 
J 


A. T | / A 7 £l 

^NIRX ± V MNIRX 


NIRXF D 


[magi 


center, J2000 


K < 15.75 




PL01 


97.96 


4.32 


0.37 


1.16 


29±5 


0.28±0.05 


8.1 ±3.3 


PL02 


98.31 


4.59 


0.94 


1.46 


32±6 


0.33±0.06 


5.1 ±2.8 


PL03 


98.38 


4.00 


0.32 


1.69 


80±9 


0.44±0.05 


6.8 ±3.2 


PL04 


98.53 


4.42 


1.10 


1.85 


89±9 


0.24±0.03 


9.7 ±5.2 


PL05 


98.63 


4.32 


0.86 


1.31 


57±8 


0.18±0.03 


7.2 ±3.0 


PL06 


98.66 


4.21 


0.73 


0.75 


13±4 


0.52±0.16 


10.4±4.0 


PL07 


98.88 


3.98 


0.38 


0.88 


22±5 


0.61±0.14 


8.9 ±3.2 


REFL08 


98.56 


4.32 


0.99 


1.30 


49±7 


0.33±0.05 


11.8±4.8 


REFL09 


98.78 


3.69 


0.74 


1.49 


65±8 


0.76±0.09 


9.8 ±5.5 


REFL10 


97.78 


5.27 


1.19 


1.15 


15±4 


0.32±0.09 


3.2 ±2.1 


NGC 2237 


97.59 


4.93 


1.94 


1.91 


36±6 


0.15±0.03 


3.1 ±1.7 


NGC 2244 


97.95 


4.94 


1.56 


2.30 


62±8 


0.12±0.02 


1.4 ±0.7 



a Number of NIRX stars with 10th Nearest Neighbor densities above 0.2 stars-arcmin 2 . 

b NIRX fraction with respect to total number of stars with K < 15.75 inside 0.2 
stars-arcmin -2 contour. 

c Average Extinction towards cluster line of sight (estimated from background source 
colors) . 
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Table 2. Rosette Complex: Sources with No detection in J, H — K > 1.5 



Brightness Limit 
[mag] 


Total 


Off Cloud 


In Cloud 
In Clusters Off Clusters 


Loose Fraction 
[%] 




[No. 


Sources] 




15.75 


112±11 


5 ±2 


86 ±9 


21±5 


19.3±4.6 


16.25 


173±13 


7 ±3 


131±11 


35±6 


20.7±3.6 


16.75 


264±16 


15±4 


194±14 


53±7 


20.7±2.7 


17.25 


330±18 


30±5 


224±15 


72±9 


23.0±2.9 



